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SLIDE STATISTICS AND FINANCIAL RETURNS 
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Abstract. We present a new approach to financial returns based on an infi¬ 
nite family of statistics called slide statistics that we introduce. The evidence 
these statistics provide suggests that certain distributions such as the stable 
distributions are not good models for the financial returns from various secu¬ 
rities or indexes like the S&P 500 and the Dow Jones. Formally, we associate 
with any finite subset of a metric space an infinite sequence of scale invariant 
numbers pi, p 2 ,... derived from a variant of differential entropy called the ge¬ 
nial entropy. We give explicit formulas for pi and p 2 that are easily evaluated 
by a computer and make this theory particularly suitable for applications. As 
statistics for point processes, these numbers often appear to converge in simu¬ 
lations and we give examples where 1/pi converges to the Hausdorff dimension 
and we prove that pi > 0. For a uniform random variable X on [a, b] n , the ev¬ 
idence from simulations suggests that pi(X) = 1/n and P 2 {X) = — n 2 / (6n 2 ) 
which yields new tests for spatial randomness. The slide statistics describe 
continuous random variables in an entirely new way. For example, if Z is any 
normal variable then simulations suggest that pi(Z) = 4/tt and P 2 (Z) = — 1 
which provides new goodness of fit tests for normality. 


1. Introduction 

We develop new entropy based statistics pi , P 2 , ■ ■ ■ called slide statistics which 
can be computed from any sample data in a metric space. As an application, we use 
these statistics to test whether financial returns are independent observations from 
a particular distribution. For example, Figure [I] shows plots of p± against n for data 
T n consisting of n-tuples of consecutive returns regarded as a subset of the metric 
space R n with the usual metric. As can be seen, the lower curve corresponding 
to the S&P 500 is very different from the ones obtained for either the Normal or 
Laplace distributions. Any potential model for the returns of the S&P 500 must be 
able simulate the p\ curve in Figure [l] which is a new requirement that is apparently 
difficult to meet. The statistics pi contain information about the interplay between 
the dimension of a space and distributions on that space and provide a new way of 
describing continuous random variables and stochastic processes in general. The pt 
will take some time to define but we will provide explicit formulas for slide statistics 
pi and p 2 that are easily computed for any sample in a metric space. 

Our applications will focus on financial returns but many other important pro¬ 
cesses in science and mathematics also yield data in the form of a set of points in 
a metric space that must be quantified and interpreted. Fields like fractal analysis 
have developed in order to obtain dimensional information from a wide range of 
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Figure 1. The lower graph is the plot of p\(T n ) against n for the 
S&P 500 for the 5000 trading days ending December 31, 2014. The 
middle cluster of five graphs are the corresponding plots obtained 
from five different simulated sequences of returns with the Laplace 
distribution. The upper cluster of five graphs are the corresponding 
plots obtained from five different simulated sequences of normally 
distributed returns. 


real world data using quantities like the Hausdorff, information and correlation di¬ 
mensions mm- One of the problems with these traditional measures of dimension 
is that they tend to concentrate on the geometric properties of a set while ignoring 
the statistical origins of the data. We will give examples where pi(U) appears to 
converge to the reciprocal of the Hausdorff dimension when U is taken to be a a 
larger and larger random sample from a fractal. 

It is not clear in general how to assign a dimension to a point process but our 
theory and simulations suggest a purely statistical way to identify a class of them 
for which it does make sense to assign a dimension. Specifically, we think the 
dimension of a point process P should be defined to be d provided pi{P) = 1/d 
and p n {P) = (—1 ) n+1 (n — l)!(n — 1 )C( n )/d n for n > 1 in which case we say the 
process is tangible with dimension d. In the case of a tangible process, the values of 
p n {P) are entirely determined by the number d but this relationship does not hold in 
general. A possible example of a tangible process is the random generation of points 
(xi, X 2 , • • ■, Xd) in [a, 6] d , where the Xi are chosen independently and uniformly at 
random from [a, b], for which we provide evidence that pi(P) = 1/d and P 2 (P) = 
— 7t 2 /(6g? 2 ). In particular, pi and p 2 are new spatial statistics for testing hypotheses 
concerning the distribution of points in [a,b] d . 

Although we will sometimes use fractals as examples, these statistics can be 
applied to any random variable with values in a metric space. In particular, the 
slide statistics appear to describe continuous random variables X in an entirely 
new way that has nothing to do with mean and standard deviation since p n (aX + 
b) = p n {X) for a / 0 . For example, the first two slide statistics for any normal 
variable Z appear to converge to Pi{Z) = 4/7r and P 2 {Z) = —1 which provide the 
basis for a new goodness of fit test for normality. For any exponential random 
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variable W, simulation gives that pi(W) ~ 1.4590 and p 2 {W) ~ —0.7333 and we 
conjecture from these and other simulations that Pi(X) > 1 for any continuous 
real-valued random variable X. If we think of l/p\(X) as the dimension of the 
random variable X , then this last conjecture says the uniform distribution on [a, b] 
has the maximum dimension of any continuous real-valued random variable. The 
second slide statistic p 2 is negative for the examples given so far but p 2 is positive for 
the Cauchy distribution and is often found to be positive for the log daily returns 
of the Standard and Poor 500 Index which is yet another demonstration of the 
non-nornrality of these returns. 

The construction of Pi(U) requires several ideas working in concert and can be 
summarized as follows. In Section [2j we introduce a variant of the differential 
entropy called the genial entropy which is the starting point for defining the slide 
statistics. Unlike the differential entropy, the genial entropy is scale invariant and 
in Section [3] we prove it is never negative which will give a new lower bound for 
the differential entropy. Given a finite set of distinct points U in a metric space, 
we find the distance from each point to its nearest neighbour and arrange these 
distances in descending order so d\ > c ?2 > • • • > d n > 0. Let f(x) be the function 
on [0,1) whose value on [(z — 1 )/n,i/n) is d,;. Let Aft) be the area under (/(.t))* 
and let aft) denote the genial entropy of the density which turns out to be 

0 at t = 0. We then define p n (U ) to be the nth derivative from the right at 0 of 
the function aft) which is developed in Section [1] 

In Section |4l we derive an explicit formula for p\(U) and state a conjectured 
formula for p2(U). The fact that pi and p2 are easily evaluated using a computer 
makes them particularly suitable for practical applications. The results obtained 
from the simulation of p\ and p 2 in a variety of contexts are summarized in Section[5] 
and demonstrate the interplay between dimension and distribution that is captured 
by these statistics. As applications in Section [6j we consider a variety of possible 
distributions as models for financial returns and use the slide statistics to illustrate 
how poorly these models fit the empirical data. 


2. The genial Entropy 

Our starting point for the development of the slide statistics is the differential 
entropy — f fln(f)dx of a density / which is well known in statistics and informa¬ 
tion theory [7]. Our focus will be on a variant of the differential entropy called the 
genial entropy or g-entropy which will be described in Definition [5] In the simplest 
case of a probability density / that has an inverse function such as —ln(x) on (0,1], 
the genial entropy is just the sum of the differential entropies of / and / -1 . The 
next proposition shows how this sum can be written in a form that makes sense for 
densities that may not have inverses. 

Proposition 1. Suppose f is a continuous function on [0,6] with the properties 
that the derivative of f exists and is negative on (0,6), /(6) = 0 and fdx = 1. 
Also assume the differential entropies of f and / -1 both exist. Then the sum 
of their differential entropies is given by — fln(f)dx — f~ 1 ln(f~ 1 )dy = 

-1 “ Jo fln(xf)dx. 
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Proof. Substitute y = f(x) into the integral — / 0 "^ 0> f~ 1 ln(f~ 1 )dy to get — / fe ° xln(x) f (x)dx 
which equals / 0 & xln(x)f'(x)dx. After integrating by parts this last integral becomes 

— 1 — Jq fln(x)dx and the result follows. □ 

The genial entropy will only be defined for densities of the following special form. 

Definition 1. A corner density is a function f : I —¥ [0, oo) where I U {0} is a 
connected interval contained in [0, oo), / is monotone decreasing and f T fdx = 1. 

Here is our definition of genial entropy which is motivated by the conclusion of 
Proposition [l] 

Definition 2. Let f •. I —¥ [0,oo) be a corner density. The genial entropy or g- 
entropy of f is defined by G(f) = — 1 — fj fln(xf)dx when this integral exists, with 
the usual convention that 01n(0) = 0. 

If / satisfies the conditions of Proposition^ then G(f) = G(/ _1 ) so in particular 

— ln(:r) and e~ x must have the same genial entropy which happens to be Euler’s 
constant as shown in Table [l] In Section [3j we will prove the genial entropy is 
always nonnegative. 


The genial entropy of some corner densities. 


Density 

Domain 

Genial Entropy 

l/b 

[0,6] 

0 

— ln(a;) 

(0,1] 

7 

e~ x 

[0, oo) 

7 

a/(x 1 ~ a ) for a £ (0,1) 

(0,1] 

— ln(a) 

2e~ x */y/n 

[0, oo) 

(—1 + 7 + ln(7r))/2 

2/(tt(1 + x 1 )) 

[0, oo) 

-1 + ln(2) + ln(7r) 


Table 1. 


Unlike the differential entropy, the genial entropy is invariant under changes of 
scale. 

Proposition 2. Let f : I —> [0, oo) be any corner density. Then for any A > 0 the 
function h : XI —► [0, oo) defined by h(z) = y/(f) for z £ XI is a corner density 
with the same genial entropy as f. 

Proof. G(h) = —1 — J xj hln(zh)dz = —1 — f XJ jf(j)ln(zjf(f))dz. After substi¬ 
tuting z = Xx , this last integral becomes —1 — J f fln(xf)dx = G(f). □ 

In Section [5j we will consider data consisting of a finite set of distinct points 
in a metric space. The first steps in computing the slide statistics will be to find 
the distances di from each point to its nearest neighbour, arrange them in decreas¬ 
ing order and use them to construct the corner density fu* discussed in the next 
definition and theorem. 

Definition 3. Let D be any sequence g?i, d 2 , ■ ■ ■, d n with di > c ?2 > • • • > d n > 0 
and let D* be the sequence d\/p, dif\x, ■ ■ ■, d n /p, where y. is the mean of the di. 

(1) Define fn ■ [0,1) —> [0, oo) to have the value di on the interval [(i — 
1 )/n,i/n). In particular, fjj* is the corner density whose value on [(?' — 
l)/n, i/n) is di/y. 
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(2) Define Lb : [0, oo) -A [0,1] by Lb{v) = m/n where m is the number of 
elements of D that are less than or equal to y. In other words , Lb is the 
restriction of the empirical cumulative distribution function for the data in 
D to the interval [0, oo). 

The next theorem shows that the genial entropy of /d* can be computed using 
the empirical cumulative distribution function. 

Theorem 1. Let D be any sequence g?i, c^, ..., d n with d\ > c ?2 > • • • > d n > 0. 
Then fB* and 1 — Lb* are corner densities with the same genial entropy. 

Proof. The sequence D may contain repetitions, so assume that the discontinuities 
of fB• in (0,1) occur at with t\ < t '2 ■ ■ ■ < i m _i and let to = 0 and t m = 1. 
Suppose the value of fB* on [tj_i, tf) is for i = 1,2,..., m and let e m +i = 0 so 
ei > e 2 > • • • > e m > e m +i = 0. We can then describe 1 — Lb* as the function 
that takes the value tj on [e*+i, ef) for i = 1,2,..., to and the value 0 on [ei, oo). 
Since fg fB*{x)dx = XaLi^/V)(*/ n — (* ~ l)/n) = 1> we must also have that 
Xu= i e i{U - U- 1 ) = 1- Rearranging this last sum gives \L{ e i - e i+ i) = 1 
so f 0 °° 1 - L d * {y)dy = 1 and 1 — Lb* is a corner density. To see that fB* and 
l— Lb* have the same genial entropy, we use to = 0, e m +i = 0 and the convention 
01n(0) = 0 to obtain: 

nl m - t . 

-1 - / fs*(x)ln(xf D *(x))dx = -1 - / eiln(xei)dx ) 

J 0 i J ti —i 

m 

— ^ ^ (^2— lCjl'fl'ij'i— l^z) 
i =1 
m 

(!) = y^(ej + itjtn(ej + itj)) - e l t i br(e i t I ) 

i=l 

= -1 - ( / tiln(yti)dy ) 

j=l “'e.+i 
/>00 

= -l- / (1 - L D *(y))ln(y(l - L D .(y)))dy 

Jo 

□ 

In view of Proposition [2j we could have defined fB* on any interval [0, 6) and 
obtained a function with the same genial entropy. The interval [0,1) was chosen 
in particular to insure the relationship between the genial entropies of fB* and 
1 — Lb* stated in part (2) of Theorem [I] We now show that the genial entropy is 
never negative. 


3. The genial entropy Inequality 

Our goal in this section is to prove that the genial entropy of a corner density 
can never be negative which we will use to prove that the first slide statistic p\ > 0. 
The idea is to first prove the necessary inequalities for step functions and then use 
the fact that monotone functions can be uniformly approximated by step functions. 
We begin with an inequality for step functions. 
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Lemma 1. Suppose that y\ > y 2 > ■ ■ ■ > y n . > 0 and 0 = to < t\ < t 2 < • • • < t n . 
Let A = 22?= i ViiU - U- 1 ). Then 22i=i(ti-iyiln(ti-iyi) - > -AinA 

Proof. Given vectors a, b £ R n , recall from OH that a majorizes b provided 
yh —1 O(j) > E i= i &(■;) for all k = 1... (n— 1) with equality for k = n and where the 
components of a and b have been sorted in descending order so am > a^ 2 ) > • • • > 
a,( n ) and 6(i) > 6( 2 ) > ••• > b( n y Alternatively, a majorizes b if vector b can be 
obtained from vector a by a sequence of ’’transfers” that allow us to change a vector 
a = (ai, a 2 ,..., ai ,..., aj ,..., a n ) into a vector b = (ai, a 2l ..., a* + A,..., dj — 
A,..., d n ) provided a* < dj and A < dj — Oj. 

We now show that the vector (A, 2 / 2 ^ 1 , 2 / 3^25 • • ■, 2/rAn-i) majorizes 
(yiti,y 2 t 2 ,... ,y n t n ) so the required inequality follows immediately from Kara- 
mata’s Inequality im and the convexity of x\nx on [0, 00 ). Since A = y n (t n — 

t n - i) + 2??=i 2 H{U - L-i) and 22?Zi 2 k{U ~ U- 1 ) > 2 hi {U ~ 1 ) = Vnt n -i, 

we can transfer an amount A = y n (t n — t n —±) from the first to the last entry 
of V\ = (A, y 2 t\, 2 / 3 ^ 2 , ■ • •, yntn-i) which means that V\ majorizes V 2 = (A— 

A, 2 / 2 * 1 ,2/3*2, • • • ,2/nin-i+A) = (E:^ 1 yiif i—ti-i), 2 / 2 ^ 1 ,2/3^2, • • • ,V n t n )- Now trans¬ 
fer an amount A = y n _i(t n -i — t n _ 2 ) from the first to the 2nd last entry of V 2 and 
continue similarly with the other entries to obtain the required majorization. □ 

Lemma 2. Suppose that 2/1 > 2/2 > • ■ ■ > y n > 0 and 0 = to < a = t\ < t 2 < • • • < 
= b. Let f be the function whose value is yi on [L-i,ti) and let P = f 0 f = y\d 

and Q = f. Then — / Q fc /(I + ln(xf))dx > PlnP — (P + Q)ln(P + Q) 

Proof. Let g be the constant function whose value is C > 0 on [it, v\ where 0 < 
u < v. Then — J u g(l + ln{xg))dx = ( Cu)ln(Cu ) — (Cv)ln(Cv) with the usual 
convention that 0 In 0 = 0 in the case when u = 0. Then — /(1 + ln(xf))dx = 

- Jo f{l + ln(xf))dx + / 0 a f{l + ln(xf))dx = 22?=i(yi t i-i ln (yi t i-i)-yitiln(y i t i )) + 
PlnP > —(P + Q)ln(P + Q) + PlnP , after applying Lemma[I] with A = P + Q. 

□ 

Lemma 3. Let f be a positive monotone decreasing function on [a, b] where 0 < 
a < b . Let P = af{a) and Q = fdx. Then —/^/( 1 + ln{xf))dx > PlnP — 
(P + Q)ln{P + Q). 

Proof. By [20], there exists a sequence of monotone decreasing step functions fk 
that converge uniformly to / with 0 < f(b) < fk{x) < /(a) for every k and every 
x £ [a, b\. By standard results, the sequence /*,( 1 + In(xfk) converges uniformly 
to /(I + ln(xf) on [a, 6]. Let Pk = afk{a) and Qk = ' fkdx. By Lemma 2 

-/ a fc / fc ( 1 + ln(xf k ))dx > P k lnP k - (P k + Qk)ln(P k + Qk) and the result now 
follows by taking the limit as k —> oo of both sides of this inequality. □ 

Lemma 4. Let f be a positive monotone decreasing function on (0, 6] for which 
Q = Jq fdx is finite. Then — / Q b /(I + ln(xf))dx > —QlnQ. 

Proof. Suppose a is between 0 and b and let P a = af(a ), Q a = fdx. Since Q 
is finite, lim a ^ 0 +P a = 0. By Lemma 0 ~ fa ft 1 + ln{xf))dx > PalnPa - (P a + 
Q a )ln(P a + Qa) and the result follows by taking the limit as a —> 0 + on both sides 
of this inequality. □ 


SLIDE STATISTICS AND FINANCIAL RETURNS 


7 


Lemma 5. Let f be a positive monotone decreasing function on (0,oo) for which 
Q = fdx is finite. Then — / 0 °° /(I + ln(xf))dx > — QlnQ. 

Proof. Let Q t = f* fdx for t > 0,. By Lemma|4j — f* f(l+ln(xf))dx > — Q t ln(Q t ) 
and the result follows by taking the limit as t —► oo on both sides of this inequality. 

□ 

Theorem 2. (The Genial Entropy Inequality) For any corner density f : I —t 
[0,oo), -1 - fj f(ln(xf)dx > 0. 

Proof. Follows from Lemma [5] by taking Q = 1. □ 

If the density function of a random variable happens to be a corner density, then 
the Genial Entropy Inequality gives a new lower bound for the differential entropy. 

Theorem 3. Let X be a random variable whose pdf is a corner density f : I —> 
[0, oo) and let h(X ) = — J fln(f)dx be the differential entropy of X. Then h(X) > 
l + _E(lnX). 

Proof. Follows from Theorem [2] □ 

4. The Slide Derivatives 

With each corner density, we now associate a function called a slide function that 
describes how the genial entropy changes as the density is deformed to a constant 
function. In Section [5j the slide numbers will be defined as the derivatives of a 
particular slide functions at 0. 

Definition 4. Suppose f : (0, b) -A (0, oo) is monotone decreasing and Aft) = 
fg(f(x)Ydx.Then <Xf(t) = = —1 — fg ^A(t) ^ n ( x ^A(t) x w be called 

the slide function of f and its domain will be the set of all t > 0 at which A and 
tjf both exist. 

By the Genial Entropy Inequality in Theorem [2j we always have agft) > 0. Note 
also that ct/( 0) = 0 and if Jq f = 1 then er/(l) = G(f). The next theorem says 
that the function af(t) is invariant under changes of scale. It follows from a simple 
change of variables argument similar to the proof of Proposition [2] 

Theorem 4. Suppose f : (0,6) -A (0, oo) is monotone decreasing and g : (0, A6) —► 
(0, oo) is defined by g(z) = /3/(f) for some A > 0 and any j3 > 0. Then <jf = a g . 

We now show that under mild conditions there must be an s > 0 for which the 
interval [0, s] is contained in the domain of Of and furthermore that Of must be 
continuous from the right at 0. 

Lemma 6. Suppose f : (0,6) -A (0,oo) is monotone decreasing with f^(f(x)) s dx < 
oo for some s > 0. Then <jf is defined on [0, s/2] and continuous from the right at 

0. 

Proof. We can assume 6 = 1 by Proposition^ Then for t £ [0, s], the function Aft) 
in DefinitionEjis finite since Aft) = J^(f(x)) t dx < fg l+(f(x)) s dx = 1+a1(s) < oo. 
To see that Um t ^o+ A(t) = 1, choose c and d with 0 < c < d < 1 and consider 
I A{t)-1\_ = | Jq (f{x)Y — ldx\ < fg \(f(x)Y-l\dx < f Q c (l + f(x) s )dx + f c d \(f(x)Y~ 
l\dx + J. (1 + f(x) s )dx. In this last sum, the first and third terms are independent 
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of t and can be made as small as desired by choosing c and d appropriately. Since 
| (f(x)y — 1 | converges uniformly to 0 on [c, d] as t —> 0 + , the results follow. 

We now show that 07 (f) is finite for t £ [0, s/2] and continuous from the right 
at 0 as follows: 

0 < 07 (f) 


= -1 - 


= -1 - 


m 


r 1 

(/ {x)) t ln{x)dx — 


m 


(f(x)) t ln((f(x)) t )dx 


A(t) 


(f(x)) t ln(A(t))dx 


= + A W ln ( A W) - J o ((fix))* - 1 )ln(x)dx 

- [ (f(x)) t ln{{f(x)) t )dx) 

Jo 

< ^yd 1 - A (t) + A (t)H A (t)\ + J o !(/(*))*-!II ln(a:)|da: 


+ [ \(f(x)YH{f(x)) t )\dx) 

Jo 

It remains to show that each of the three terms in this last sum is finite for 
t € [0, s/2] and goes to 0 as t -A- 0 + . For the first term, clearly 
lim t _ >0 + (l — Aft) + A(t ) ln(H(f)) = 0 . 

For the second term, f Q \(f(x)Y — l||ln(a;)|dx converges since the integrand 
is the product of square integrable functions. To see that lim t _ >0 + | (f(x))* — 
l||ln(a:)|dx = 0 , choose c and d with 0 < c < d < 1 and consider the inequality 
/o I (f(x)Y - 1 || ln(x)\dx < f Q c ( 1 + f(x) s / 2 )\ ln(a;)|da: + \(f(x))* - 1 || ln(a;)|dx + 

J, (1 + f(x) s / 2 ) | ln(a;)|dx. Now follow the argument given above for |A(t) — 1|. 

For the third term, use the inequality z — 1 < zlnz < z(z — 1) for 2 > 0 to 
get f* \(f(x)) t ln((f(x)) t )\dx < fg |(/(x)) t — 1 ] max(l, {f{x)) t )dx which converges 
since the last integrand is the product of square integrable functions. Now show 
lim t ^ 0 + Jq \(f(x)Y — 1 | max(l, ( f{x)) t )dx = 0 as before. 

□ 


The information we wish to extract from a corner density is captured by the 
derivatives of its slide function at 0 which we now describe. 

Definition 5. Suppose f : (0, b ) -A (0, 00 ) is monotone decreasing with j^(f(x)) s dx < 
00 for some s > 0. Then the n’th slide derivative of f is defined by ip n (f) = d j£/ (0) 
where all derivatives are taken from the right. If all of these derivatives exist, then 
the slide series of f is defined to be —“r^f"- 

Here are some elementary properties of ipn(f)- 

Theorem 5. Suppose f : (0, b) -A (0, 00 ) is monotone decreasing with J^(f(x)) s dx < 
00 for some s > 0 . 

( 1 ) Ifipi(f) exists then ipi(f) > 0- 
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( 2 ) Ifipn(f) exists then so does if n (f r ) for r > 0 and if n (f r ) = r n ip n {f). 

(3) If f is a constant function, then ip n {f) = 0 for all n. 

Proof. (1) The slide function tjf is non-negative on its domain and er/(0) = 0 so 
the first derivative must be nonnegative. 

(2) For t sufficiently small and nonnegative, we have <Jfr(t) = <Jf(rt ) and the 
result follows from the chain rule. 

(3) If / is a constant function, then <r/(t) = 0 for all t > 0. 

□ 


Here is an example of a slide derivative calculation that we will connect with the 
uniform distribution on [0,1] in Section [5j 

Theorem 6. Let /( x) = — ln(a;) for x £ (0,1). Then for t > 0, cr/(f) = — 1 + t — 
t'L(f) + log(r(l + 1)) and the slide derivatives of f are given by tpiif) = 1 an d by 
i>n(f) = (-l) n+1 (™ - l)!(n - l)C(n) for n > 1 . 


Proof. Let g t (x) = = ( r ^+t) so Jo 9 t(x)dx = 1 . Then 


( 2 ) 


= G(g t (x)) 

= -1 - / gt ln(xg t )dx 
Jo 


g t ln(x)dx - g t ln(g t )dx 
Jo 

= -1 - (-1 - t) - (1 - ln(T(l + t)) + tV(t)) 
= — 1 + t + ln(T(l + f)) — 



The result now follows by differentiating — 1 + t + ln(T(l + t)) — f’F(f) from the 
right at t = 0 . □ 


Corollary 1 . Let f r {x) = (— ln(a;)) r for x £ (0,1) and r > 0. Then V’i(/r) = r 
and ifnifr) = (— l) n+1 (n — l)!(n — l)C(n)r” for n > 1. The slide series for f r is 
given by rt + £“ 2 

Proof. Immediate from Theorem [5] and Theorem [ 6 j 

□ 


Example 1. ipi(f) = oo for the function f{x) = exp(— 1/(1 — a;) 2 ) on [0,1). 

In order to apply the slide derivatives to samples in a metric space, we first find 
the distance di from the Fth point to its nearest neighbour. These distances can be 
used to construct the function fo in Definition [3] and we can then compute the first 
slide derivative ifiifo) using the explicit formula provided in the next theorem. 

Theorem 7. Suppose D = {di, d 2j ■ ■ •, d n } where we assume d\ > d 2 > • • • > 
d n > 0 and let fo be the function on [ 0 , 1 ) whose value on the interval i ) is 
di as in Definition [3| Then the first slide derivative of fo is given by V'i(/d) = 

n T,i=2 * ln (*) ^ ElL'i 1 M^). 

Proof. By Definition [5j we have to calculate the right hand derivative of the slide 
function at t = 0 so we first find an expression for <Jf D (t). Let gfix) = 
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w h ere A(f) = Jq (foixjydx so gt{x) takes the value a,;(t) on [^4, where 


E"= i ' 


o-/(0 = G (9t(x)) 


= -l- g t (x)ln(xg t {x))dx 
Jo 

n p i- 

= -1 - 5Z i _ 1 


= -l- ((^)ln(^) - (^) +y ((^i)ln(^l) - ( (i 1)Q?; )ln( (i 1)0i ) - *)) 
vv n y v n ' v n y ' vv n y v n ' v n ' v n y n 

i =2 


-'flu , /« 1 \ , /iO-i\ V - ^ < (* ~ l) a i \ l / (* — l)<2iN 


~(^> ta <^>-E(=)M=)+E( 


n n 


To find the derivative of <Tf(t), we now use the facts that ln(d/c) — 


£ £!W M«y aiK| AlA = f 1 + i"«)S • 




( i _iK) ln (( i _iK)) ( _o 


= -( 1 + 1 ”(^44))(^)(( = A) 1 ”( <i i)- ;E ta (<ij)) 


- £ ((!+1»(^))(^)(( V> ‘"W - s E‘“W») 

i=2 j#* 

+ E ((1 + (1^M))(^ )( (^) ln(*) - i Eln( dj ))) 


At t = 0, each a* is equal to 1 so this derivative becomes 
^(0) = -(l + ln( 1 ))( 1 )((-— 1 ) ln(di) - 1 ^Sn(og) 

at v Ti Tl Tl 71 


- E ft 1 +1“ (b) 0« V> '“W -; £ ^w))) 

i=2 j^i 

n ; _ 1 : _ 1 _ 1 -1 n 
+ ((! + ln ( — ))(—)((——) ln K) - - 5Z ln K'))) 

i=2 ' ' 

= - (! + ln (^)) (^) ((^) ln ( dl ) - ^ E ln (^)) 


+ 51 ((-! + (* - l)ln(i- 1) -*ln(i) +ln(n))(t)((^—i)ln(di) - t ^ln(dj))) 
*—2 

= Pi + P 2 + P 3 
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The Pi terms are defined and calculated as follows. The P\ term consists of the 
parts of this last expression that involve the isolated Is. 

1 1 1 ^ 

Pl = "Wb) ((^H ln(dl) ~ n E ln ^)) 

1 

n 1 _ 1 1 n 

+e ((-!) y ((V) w ^ e M^))) 

*= 2 

= 0 


P 2 is the sum of all of the terms containing Inn. 


p 2 = (^r) (( n " x ) ln ( d l) - E ln ^)) 


+ (-y e (((«- ^ in (^) -E in (^))) 

i=2 

+ ( ”? n )(( n - !) ln ( d «) - E ln ( rf j)) 

= ( E nn )((»^l) ln(d„)-^ ln (^)) 

= ^)v ln( y 

n f-f d n 

l—l 


P 3 is what remains after Pi and P 2 are subtracted from 0). 

n _ 1 ^ n 

p 3 = E (((* -!) ln (* - x ) - * ln (*)) (-) ((——) H d i ) - - E ln ( d f))) 


i —2 




j #*+1 


+ (^2“)(( n “ !) ln ( d «) - E ln ( rf j)) 

1 n—1 n n 

= “2 E* lni (( ~ ( n ~ 1 ) ln ( rf i) + E ln ( d J')) + i) ^ ln(dj))) 

<=2 

^ n—1 

= — ^ f lni( - nln(d») + nln(d i+ i)) 

n i=2 

1 n-1 u 

= -E^WM-x 1 ) 


i=2 


□ 


The following formula for the second slide derivative if >2 (/) is motivated by cal¬ 
culations for small values of n. In the next section, we will see that the results from 
simulations based on this formula agree with what we would expect from theoretical 
considerations. 
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Conjecture 1. Suppose d\ > d 2 > • • • > d n > 0 and let fr> be the function on 
[0,1) whose value on the interval [^=1, -) is dj. Let S/ = ]C" =1 l°g(di), &2 = 
X)"=i l°g(^i) 2 an ^ = Yli=i l° 9 (di/dn) 2 ■ Then the second slide derivative of fjj 
is given by 

n— 1 

MId) = -(5Z (*log(«) \og(d i+ i/di){2Si - nlog(djdj + i))) 
i=1 

+ log(n)(2(5i - nlog(d„)) 2 - nS 3 ) + nS 2 - S{)/n 2 . 

The next section define the slide numbers and illustrates their application to 
some standard point processes. 

5. The Slide Numbers 

Sample data often consists of a set of distinct points in a metric space. With the 
help of Definition [5j we can now associate with each of these samples an infinite 
family of new statistics called the slide numbers. 

Definition 6. Let M be a metric space and let U = {zii,zi 2 ,... ,ztfc} be a set of k 
distinct points in M. For each i = 1,..., k, let di be the distance from Ui to its 
nearest neighbour in U. Define a sequence D by ordering the di in descending order 
as dji] > dp] > • • • > d[k] > 0. ^4s in Definition [3j let fu be the function on [0,1) 
whose value on the interval is duy Define the n’th slide number of U by 

p n {U) = 'ipnifo) and define the slide series of U to be ^pp-t n . 

Values of pi(U) for various random variables are shown in Table [2] which shows 
the connection between l/pi(U) and the Hausdorff dimension of [0, l] m , the Cantor 
set and the Sierpinski triangle. For certain point processes P, the statistics p n {U) 
appear to converge as the sample size U gets large. For example, for any normal 
random variable Z 1 the quantity pi(U) appears to converge to 4/-7T in which case it 
makes sense to write pi{Z) = 4/7r. More generally we define p n {P) for an arbitrary 
point process P as follows. 

Definition 7. Let M be a metric space and let U = {ui,U 2 , ■ • • } be a sequence of 
distinct points in M generated by some point process P and letUk = {zti, «2> • • ■, Wfc}- 
If p n (Uk) converges in probability as k -A oo, then p n (P) is defined to be the value 
of this limit. If all of the limits p n (P ) exist, then we define the slide series of the 
process P to be t n . In the case where U is a sample of a random variable 

X, we will use the notation p n (X) instead of p n (P). 

The following proposition follows immediately from the definitions and shows 
that adjusting the mean or standard deviation of a random variable X has no 
effect on p n (X). 

Proposition 3. If X is a random variable, then p n [aX + b) = p n {X) for all real 
number a and b with a/0. 

Some evidence for the convergence of the slide statistics p n {U) is given in Table[2] 
and Table [3] In these tables, the points in the Cantor set were generated using 
fi where the cq are either 0 or 2 with probability 1/2. Points in the Sierpinski 
triangle were generated using the Chaos Game pQ. The generation of all random 
numbers used in these simulations was based on the Mersenne Twister. 
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Simulated Values of pi(U) for Various Random Variables 


Density 

Mi 

Ol 

1 /Mi 

conjectured 

uniform on [a, b] 

1.0003 

0.0111 

0.9997 

i 

normal 

1.2664 

0.129 

0.7896 

tt/4 ~ 0.785 

exponential 

1.4590 

0.0141 

0.6854 

? 

l/(2y/x) on [ 0 , 1 ] 

1.2817 

0.0132 

0.7802 

? 

uniform on [ 0 , l ] 2 

0.5023 

0.0056 

1.9908 

2 

uniform on [ 0 , l ] 3 

0.3416 

0.0037 

2.9274 

3 

uniform on [ 0 , l ] 4 

0.2642 

0.0029 

3.7850 

4 

bivariate normal 

0.7264 

0.0073 

1.3766 

? 

Cantor 

1.6014 

0.0170 

0.6244 

ln(2)/ln(3) « 0.631 

Sierpinski 

0.6344 

0.0067 

1.57624 

ln(3)/ln(2) « 1.5849 


Table 2. For each density, 1000 samples of size 10000 were 
generated and the value of P\{U) was computed for each sample 
U. The mean pi and standard deviation < 7 i of these 1000 values 
for p\{U) are shown. The value given for 1/pi is the conjectured 
limiting value of l/pi as the sample size approaches infinity. 


Simulated Values of P 2 {U) for Various Random Variables 


Density 

M 2 

<J2 

P 2 conjectured 

uniform on [a, b] 

-1.6461 

.0732 

-tt 2 /6 ~ -1.6449 

normal 

-1.0273 

0.0860 

-1 

exponential 

-0.7333 

0.0920 

? 

l/( 2 y / x) on [ 0 , 1 ] 

-2.5792 

0.1085 

? 

uniform on [ 0 , l ] 2 

-0.4096 

0.0186 

-(tt 2 /6)(1/2 ) 2 « -.4112 

uniform on [ 0 , l] a 

-0.1825 

0.0083 

-(tt'76)(1/3 ) 2 « -0.1827 

uniform on [ 0 , l ] 4 

-0.1038 

0.0049 

-(tt 2 /6)(1/4 ) 2 « 0.1028 

bivariate normal 

-0.2004 

0.0233 

? 

Cantor 

-4.1464 

0.1933 

(_l)^ + i(2-l)!(2-l)C(2) _ , i nr, 

(ln(2)/ln(3)) 2 ~ T.IOZ 

Sierpinski 

-0.6549 

0.0295 

(_l)^+i(2-l)!(2-l)C(2) _ n 

(ln(3)/ln(2)) 2 ~ U ' U0 ° 


Table 3. For each density, 1000 samples of size 10000 were 
generated and the value of P 2 {U) was computed for each sample 
U using Conjecture [l] The mean p 2 and standard deviation eg of 
these 1000 values for p 2 {U ) are shown. The value given for p 2 is 
the conjectured limiting value of p 2 as the sample size approaches 
infinity. 


Consider the following thin outline for a possible argument to explain the empiri¬ 
cal results obtained for [0, l] m in Tableland Table[3] Suppose U = {zii, U 2 , ■ ■ ., Uk} 
is a very large sample of points (x±,X 2 , ■ ■ ■, x m ) in [ 0 , l] m where the ay are chosen in¬ 
dependently and uniformly at random from [ 0 , 1 ], and let u be a particular point in 
U. By (21] . the probability that a point is within r of u is approximately 1 — e~ ar 
for some a. If the sample U is large enough, then the set of nearest neighbour 
distances will be sufficiently independent mm that their empirical cumulative 
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distribution will also be approximately equal to 1 — e~ arm . If D is the ordered 
sequence of nearest neighbour distances, then Lp* ( r ) in Definition [3] will be ap¬ 
proximately 1 — e -/3r for some /3 and 1 — Ljj*(r) will be approximately e~ Sr 

JL. 

with inverse fD*(x) = ■ Now fo*{x) has the same genial entropy as 

1 — Ljj*{r) by Theorem [l] so Corollary [l] now suggests the following conjecture 
which is supported by the empirical results shown in Tableland Table [3] 

Conjecture 2. Let U = {u\,U 2 , ■ ■ ■) be a sequence of points (aq, X 2 , ■ ■ ■, x m ) in 
[0, l] m , where the Xi are chosen independently and uniformly at random from [0,1] 
and let 14 = {iq, iq, ■ ■ ■, Wfe}- Then as k —> oo, pi(£4) converges in probability 
to 1/m and p n (Uk) converges in probability to (—1 ) n+1 (n — l)!(n — 1 )((n)/m n for 
n > 1. 


There appear to be cases other than [0, l] m for which the dimension equals m 
and p n converges to (—l) n+1 (n — l)!(n — l)f(n)/m n . For example, if we take 
m = log(2)/log(3) then (—1) 2+1 (2 — 1)!(2 — l)£(2)/m 2 ~ —4.132 which is close 
to the value shown in Table [3] for the Cantor set. A similar result holds for the 
Sierpinski triangle and prompts us to make the following definition. 

Definition 8 . In the context of Definition we say that a point process P is 
tangible provided there is a number d with pi(P) = 1 /d and p n (P ) = (—1 ) n+1 (ro — 
l)!(n — 1 )(,{n)/d n for n > 1. The number d will be called the slide dimension of the 
process. If there is no such number, the process will be called intangible. 


The values for pi and p 2 given in Table [2] and Table [3] suggest that the normal 
distribution does not satisfy the conditions for tangibility in Definition [ 8 ] so cannot 
be assigned a slide dimension. In the case of a tangible process, we can recover 
the dimension by rearranging the second derivative to obtain d = ■ This 

relationship is particularly useful for spatial statistics because it provides consider¬ 
ably better estimates for the dimension of [ 0 , l] m than the values for l//q shown in 
T able [2) For a tangible p rocess P , we would like to know in general if the statistics 

( 1 ) + converge more quickly to the dimension for larger values 


of n. 

For the subset U of R generated by 20000 iterations of aq+i = Xi + cos{i) with 
Xq = 0, the value of p\{U) is approximately 0.53. But values of p\{U) less than 
1 cannot occur for continuous real-valued random variables according to the next 
conjecture which is supported by the results in Table [2] 


Conjecture 3. If X is a continuous real-valued random variable for which pi(X) 
exists, then pi(A') > 1. 


If U consists of the first 20, 000, 000 primes, then the value of p\(U) is approxi¬ 
mately 0.77235 which is interesting in view of Conjecture [3j 

The values of P 2 {X) are all negative in Table[3]but P 2 {X) appears to be positive 
when X has a Cauchy distribution which raises the following question. 


Question 1. What conditions on a random variable X determine the sign of p n (X) 
when n > 1 ? 


6. Applications of the slide statistics to financial returns. 

If the daily closing prices of an index or stock are Xi,X 2 , ■.. and the daily returns 
are given by Ui = \n{X i+ i/Xf), then a central question of financial mathematics is 
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the problem of describing the distribution of the returns U, . We now show that the 
slide statistics impose strict constraints that allow us to quickly reject many possible 
candidates for this return distribution. Our approach will be to fix a sample size r 
and use the standard metric on R n to calculate p± and P 2 for the subset of R n given 
by T n = {(Ui, U 2 ,U n ), {U 2 ,U 3 ,..., U n+1 ),(U r , U 3 ,..., U n+r -i)}. Note that by 
the same reasoning used in Proposition [3j the values of pi(T n ) are not changed if we 
replace the U by aUi + b with o / 0. In other words, Pi(T n ) is detecting information 
about the returns that has nothing to do with their mean and standard deviation. 

Figure [I] shows that the values of pi(T n ) obtained for the S&P 500 are very 
different than the values obtained for either the Normal or Laplace distributions. 
The lowest of the 11 curves shown in Figure [l] is the plot of pi(T n ) against n 
for the S&P 500 for the 5000 trading days ending December 31, 2014. The middle 
cluster of five curves are the corresponding plots obtained for five different simulated 
sequences of returns having the Laplace distribution. The upper cluster of five 
curves was obtained for the Normal distribution and we see that these curves are a 
bad fit to the p\ curve for the S&P 500. Any suggested model for the returns of the 
S&P 500 must be able to generate an approximation of this p\ curve which places a 
tight constraint on potential stochastic processes. In Figure|2j the values of p 2 {T n ) 
shown on the upper graph were obtained for the S&P 500 for the 5000 trading days 
ending December 31, 2014. For this particular time period, p 2 (T n ) was positive but 
the values obtained for samples from Normal distributions are typically negative as 
shown by the 5 graphs below the horizontal axis in Figure [2j 



n 


Figure 2. The lower graph is the plot of pi(T n ) against n for the 
500 monthly returns of the S&P 500 up to the end of 2014. The 
upper cluster of five graphs are the corresponding plots obtained 
from five different simulated sequences of normally distributed re¬ 
turns. 

The failure of the daily returns to be normal might be due to dependencies in the 
consecutive daily returns and it is possible that we might obtain a better fit using 
the monthly returns instead. The bottom curve in Figure [2] shows p\{T n ) plotted 
against n for the 500 monthly returns of the S&P 500 up to the end of 2014. The 
group of five curves at the top of Figure [2] are the corresponding plots obtained for 
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five different simulated sequences of returns having a Normal distribution. The fit 
is better than what we obtained for the daily returns but once again the normal p\ 
curves are not a match for the p\ curve for the S&P 500 which is consistent with 
the non-normality of such returns found in [5]- As another point of view on the 
non-normality of the monthly returns, Figure 3] shows the points (pa (T 3 ), p\ (X 3 )) 
plotted for 100,000 samples from a normal distribution as well as the point obtained 
for the 500 monthly returns of the S&P 500 up to the end of 2014. The possibility 
of describing these monthly returns as a mixture of two Gaussians is considered in 
[5] but we find P 2 (T 3 ) < 0 for their models but we sometimes find p‘i{T\i) > 0 for 


the S&P 500. 



Pi 


Figure 3. The small square at the upper left plots the point 
(p 2 (T 3 ), pi (T 3 )) obtained for the 500 monthly returns of the S&P 
500 ending with December of 2014. The dots are the points 
(p 2 (X 3 ), pi (X 3 )) obtained from 100,000 different simulated se¬ 
quences of normally distributed returns. 


Values of pi and P 2 for Various Stocks and Indexes 


Name 

piCt 3 ) 

P2CP3) 

Exxon 

1.44659 

-0.12179 

FTSE 

1.34232 

-0.08638 

IBM 

1.32328 

0.011604 

3M 

1.38742 

-0.07126 

NASDAQ 

1.24174 

0.03204 

Russell 2000 

1.33124 

-0.12882 

S&P 500 

1.41405 

-0.02228 


Table 4. Each of these values was calculated using the 5000 
trading days ending with June 30, 2008. 


The stable distributions are often considered mmm as possible models for 
financial returns. We now consider the family of stable distributions S(a, /3, 7 , J) 
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n 


Figure 4. The upper graph is the plot of P 2 {T n ) against n for the 
S&P 500 for the 5000 trading days ending December 31, 2014. The 
lower cluster of five graphs are the corresponding plots obtained 
from five different simulated sequences of normally distributed re¬ 
turns. 

described in [ 22 ] where a £ (0,2] is the stable parameter, j3 £ [—1,1] is the skewness 
parameter, 7 > 0 is the scale parameter and S £ (— 00 , 00 ) is the location parameter. 
By Proposition [3j the values of pi are not affected by changes in 7 and S so we will 
simply work with the family S(a, /3,1, 0) which we write as S(a, /3). Figure [5] shows 
the points (^ 2 (^ 3 ), pi(T 3 )) from Table [4] plotted together with 100,000 samples 
from stable distributions with parameters a and j3 chosen uniformly at random 
from the intervals (1,2) and (0,1) respectively. We see that the values in Table [4] 
lie outside the region corresponding to the stable distributions. Figure [ 6 ] show 
an expanded view of Figure [5] showing that the slide statistics for the financial 
returns in Table [4] are well outside the region corresponding to samples from stable 
distributions. This failure of the stable distributions to fit financial returns is 
consistent with the findings of mm- 

7. Conclusions and Future Work 

As we have seen, the statistics p\ and P 2 can be used as the basis for goodness 
of fit test for financial returns. Much more research needs to be done to better 
understand the role of the slide statistics in characterizing financial data. In par¬ 
ticular, we would like to understand the relationship between the sign of P 2 and 
the behaviour of financial markets. The simulations we have described suggest 
that pi and P 2 can be used to distinguish between probability distributions and 
are capapble of detecting dimensional information. The values of p n we obtained 
through simulation point to a larger theory of these statistics which is currently 
being developed. At present however, the p n are quite mysterious and much work 
will need to be done to understand all they are telling us about sets of points in 
metric spaces, random variables or point processes in general. 

We gave examples of point processes on the Cantor set and the Sierpinski tri¬ 
angle for which 1/pi converged to the dimension of the fractal. We would like to 
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Figure 5. The dots plotted here are the points (pi(T 3 ), p 2 (T 3 )) 
obtained for 100, 000 samples from stable distributions with pa¬ 
rameters a and /3 chosen uniformly at random from the intervals 
(1,2) and (0,1) respectively. The values in Table [4] are shown 
plotted with squares. 


0.8 -i- 

o.7 -U^ ^ . ’ ; 

0.6 

0.5 

0.4 

P2 

0.3 

0.2 

0.1 

0 . 

12 1.25 1.3 1.35 1.4 1.45 1.5 

-0.2 -I- 

Pi 

Figure 6. Expanded view of Figure [5] 

understand when this occurs in general and also the relationship between l/pi and 
the usual definitions of dimension. More generally, we would like to know when a 
point process is tangible in the se nse of Definition| 8 ) W hen a process is tangible it is 

possible to use the expression \J + as a statistic for estimating 

the dimension. We gave examples in which it converged to the dimension faster for 
n = 2 than for n = 1 and we would like to know what happens for larger values of 

n. 

In terms of calculations, we need to prove Conjecture [I] concerning the calcula¬ 
tion of P 2 (U) and we need formulas for p n (U) for n > 2. Given the complexity of 
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our conjecture for p 2 (U), the formulas for larger values of n are likely to be very 
complicated. The convergence of the p n (U) for real-valued random variables can 
sometimes be improved by using the distances between consecutive points rather 
than the distances to nearest neighbours. We would like to have a better under¬ 
standing of this situation and also to know if there are any higher dimensional 
analogues. 

In defining the slide statistics, we used the functions ( f(x)Y which can be thought 
of as a continuous deformation of f(x) at t = 1 into the constant function 1 at t = 0. 
We can achieve the same effect using the functions tf(x) + (1 — t) and calculate 
the derivatives corresponding to the slide derivatives in Section [4] It turns out that 
only the first of these statistics is interesting and while it is easier to calculate than 
Pi it often doesn’t converge. 


References 

[1] Barnsley, M., Fractals Everywhere , Academic Press, 1988. 

[2] Billingsley, P., Probability and Measure , John Wiley and Sons Ltd, 1995. 

[3] Bonetti, M., Forsberg, L., OzonofF, A. and Pagano, M., The distribution of interpoint dis¬ 
tances. Mathematical Modeling Applications in Homeland Security , HT Banks and C Castillo 
Chavez, Eds., 2003,87-106. 

[4] Belov,I., Kabainskas,A.,Sakalauskas,L., A Study of Stable Models of Stock Markets, Informa¬ 
tion Technology and Control, 2006, Vol.35, No.l 

[5] Behr, A.,Ptter, U., Alternatives to the normal model of stock returns: Gaussian mixture, 
generalised logF and generalised hyperbolic models , Annals of Finance, 2009, 5:4968 

[6] Chambers, J.M.,Mallows,C.L., Stuck,B.W., Method for Simulating Stable Random Variables, 
Journal of the American Statistical Association, 1976, Vol. 71, No. 354, pages 340-344 

[7] Cover, T.,Thomas, J. Elements of Information Theory , John Wiley and Sons, 2006. 

[8] Edgar, G., Measure, Topology, and Fractal Geometry , Springer-Verlag, 1990. 

[9] Falconer, K., Fractal Geometry: Mathematical Foundations and Applications , John Wiley 
and Sons, Ltd., 2003. 

[10] Grassberger, P. P., Procaccia, I., Measuring the strangeness of strange attractors , Physica, 
9D, 1983, 189-208. 

[11] Harte, D., Multifractals, London: Chapman and Hall, 2001. 

[12] Haas,M.,Pigorsch,C., Financial Economics, Fat-Tailed Distributions, Complex Systems in 
Finance and Econometrics, 2011, pp 308-339 

[13] Kay 11, M., Integrals Don’t Have Anything to Do with Discrete Math, Do They?, Mathematics 
Magazine 84(2);2011, 108-119 

[14] Kester, A., Asymptotic Normality of the Number of Small Sistances between Random Points 
in a Cube, Stoch. Proc. Appl. 3,1975, 45-54. 

[15] Klement, E. P., Mesiar, R., Pap, E., Quasi- and pseudo-inverses of monotone functions, and 
the construction of t-norms, Fuzzy Sets and Systems, 104, 1999, 3-13. 

[16] Lau,A., Lau,H., Wingender,J., The Distribution of Stock Returns: New Evidence against the 
Stable Model, Journal of Business and Economic Statistics, 1990, vol. 8, issue 2, pages 217-23 

[17] Marshall, A., Olkin, I., Inequalities: Theory of Majorization and Its Applications, Academic 
Press, 1980. 

[18] Shurygin, A., Using interpoint distances for pattern recognition, Pattern Recognition and 
Image Analysis,16(4),2006,726-729. 

[19] Simkowitz,M., and Beedles,W., Asymmetric Stable Distributed Security Returns, Journal of 
the American Statistical Association Vol. 75, No. 370 (Jun., 1980), pp. 306-312 

[20] Sohrab, H., Basic Real Analysis, Birkhauser Boston, 2003. 

[21] Stoyan, D., Kendall, W., Mecke, J., Stochastic Geometry and its Applications, John Wiley 
and Sons Ltd, 2008. 

[22] Weron, R., On the Chambers-Mallows-Stuck method for simulating skewed stable random 
variables, Statistics and Probability Letters 28 (1996) 165-1771 


20 


WILLIAM J. RALPH 


Mathematics Department, Brock University, St. Catharines, Ontario, Canada L2S 
3A1, Email: bralph@brocku.ca, Phone: (905) 688-5550 x3804 



